Magnetic transition and orbital degrees of freedom in vanadium spinels 
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We propose a scenario for the two phase transitions in AV2O4 (A=Zn, Mg, Cd), based on an 
effective spin-orbital model on the pyrochlore lattice. At high temperatures, spin correlations are 
strongly frustrated due to the lattice structure, and the transition at ~50 [K] is an orbital order, 
supported by Jahn- Teller lattice distortion. This orbital order introduces spatial modulation of spin 
exchange couplings depending on the bond direction. This partially releases the frustration, and 
leads to a spin order at ~40 [K]. We also study the stable spin configuration by taking account 
of third-neighbor exchange couplings and quantum fluctuations. The result is consistent with the 
experimental results. 

PACS numbers: 75.10.Jm, 75.30.Et, 75.30.Ds, 75.50.Ee 



Pyrochlore lattice is a network of corner-sharing tetra- 
hedra shown in Fig. 1, and it is a typical geometrically 
frustrated system in three dimensions. The vanadium 
spinel, ZnV204, is an insulator, and its sublattice of mag- 
netic vanadium ions constitutes a pyrochlore lattice. It 
was found that this compound reveals two phase transi- 
tions at T c i=50[K] and T c2 =40[K]. 1 An X-ray diffraction 
experiment showed that the transition at T& is a struc- 
tural transition from the high-temperature cubic phase 
to the low-temperature tetragonal phase with the lat- 
tice constants a = b > c. The neutron experiment at 
T = 4.2[K] showed the presence of the antiferromagnetic 
long-range order plotted in Fig. l£ and the Li-substitute 
material showed an anomaly in its NMR signal at T&m 
These indicate that the lower-temperature transition in 
ZnV204 is a paramagnetic-to-antiferromagnetic transi- 
tion. 

Theoretical studies have predicted unusual properties 
for spin systems on the pyrochlore lattice. In particular, 
antiferromagnetic classical spin systems with only near- 
est neighbor interactions are believed to have no mag- 
netic order at any temperature^^ It is also believed 
that the quantum spin systems have a spin-singlet ground 
state and a finite energy gap to spin-triplet excitations 
with thermodynamic number of singlet states inside the 
singlet-to-triplet gap. Several symmetry breakings are 
theoretically predicted within the spin singlet subspace, 
e.g., dimer/tetramer order, but without magnetic long 
range order £2*2*i2i 

Therefore, we encounter a difficulty in explaining these 
phase transitions in Z11V2O4, if the material is consid- 
ered as a pure spin S — 1 system on the pyrochlore lat- 
tice, and other degrees of freedom are necessary to take 
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FIG. 1: Cubic unit cell of the pyrochlore lattice and spin or- 
der determined by neutron experiments at low temperatures. 



into account. Yamashita and Ueda studied this problem 
based on a valence-bond-solid (VBS) approach and ex- 
amined the effects of Jahn- Teller distortion^ They pro- 
posed that the transition at T c \ is due to the Jahn- Teller 
effect which lifts the degeneracy of the spin-singlet local 
ground states at each tetrahedron unit of the pyrochlore 
lattice. This idea was also applied to classical spin sys- 
tems, and the effects of the coupling to the lattice dis- 
tortion were investigated using the point-group argument 
and the Landau theory— These scenarios are quite ap- 
pealing, but some difficulty still seems to remain. 

The problem is that it is difficult to explain the mag- 
netic transition at the low transition temperature T C 2 by 
Yamashita-Ueda type scenarios based on a quantum-spin 
picture. Like other theoretical works, they started from 
the assumption that there exists a finite energy gap be- 
tween the spin-singlet ground state and spin-triplet ex- 
citations, and constructed a low-energy effective theory 
to describe a phase transition within the spin-singlet sec- 
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tor. High-energy excitations with total spin 5 ^ are 
already traced out from the theoretical framework, and 
there remain no degrees of freedom describing the low- 
temperature magnetic transition. 

Scenarios based on a classical spin picture do not have 
this problem, but it has another difficulty to explain the 
following generic difference between the vanadium and 
chromium spinels. The presence of these two transitions 
is common to other vanadium spinels, MgA^O^i 3 - and 
CdV204ji^ with the same valence V 3+ . On the other 
hand, the chromium spinels, ZnCr 2 04, CdCr 2 04, and 
MgCr204 show only one transition at 12.5[K], 7.8[K], 
and f 2.5[K], respectivelyii£*i& The difference between the 
vanadium spinels and the chromium spinels is generic 
and independent of divalent A-site cations. It is also 
unclear whether classical approximations are justified at 
low temperatures for the system with the second smallest 
spin, 5 = 1. 

In the present study, we will explore another scenario 
to explain the two transitions with taking account of the 
orbital degrees of freedom, which exist in the vanadium 
spinels but not in the chromium spinels. The essential 
difference between the vanadium and chromium spinels 
is the number of electrons in magnetic ions: V 3+ ions 
has d 2 configuration and Cr 3+ has d 3 configuration. Due 
to the cubic crystal field, the d-electron orbitals are split 
into the high-energy e g and low-energy t 2g multiplets. 
Since in both the spinels the high spin state is realized 
due to large intra-atomic Coulomb interactions, two of 
the 3-fold degenerate t 2g orbitals are occupied by elec- 
trons in the vanadium case, whereas all three orbitals 
are occupied in the chromium case. Therefore, each V 3+ 
ion has 3-fold orbital degeneracy in addition to the triplet 
spin state 5 = 1. Rigorously speaking, the crystal field 
has a small trigonal component, which may result in a 
further splitting of the t2 g multiplet. However, this split- 
ting is compensated by the covalency difference of the 
three states, and the net splitting will be small and we 
neglect this. 

The realistic model Hamiltonian for this system reads 
with the standard notation, 

H = S [ ta ^ - r Mcca C 3fl° + H ' C - 

a/3 a 

+ ^Y Y Y, Ua f } > a 'P' C ic«r C 0T C i0''r c ia>a, (1) 
i a{3,a'f3' err 

where i,j are site index, a, t are spin index, and a, /3 = 
1 (d yz ), 2 (d zx ), 3 (d xy ) are orbital index. For the 
Coulomb interactions, we use the standard parameter- 
ization, U a p tOL ipi = VS aa '5fjp' + J{5 a f]>5p a > + 5 a p8 a ipi), 

and U = V + 2 J is the Hubbard interaction in the same 
orbital. These values were determined from the spec- 
troscopy data for the vanadium perovskites as, U ~6[eV] 
and J ~0.7[eV]fii and it is reasonable to assume a small 
value for the ratio of the two in the vanadium spinels, 
r\ = J/U ~ 0.11. The hopping integrals t a p(ri — Tj) 
were also determined by a parameter fitting of the first- 



principle band calculation data^ and it was found that 
the one for the nearest-neighbor cr-bond is much larger 
than the others. We therefore consider only this kind of 
hoppings, t a = — 0.32[cV], and set the others zero. This 
is the hopping in the case where one of four lobes of each 
tig orbital is pointing towards the other end of the bond. 
In this case, the number of electrons in each orbital is 
conserved through the hopping processes, simplifying the 
following calculations. We neglect the relativistic spin- 
orbit coupling, which is much smaller than the energy 
scales in Eq. QJ. 

Since ZnV204 is an insulator, we employ the strong 
coupling approach U 3> t a , and include the hopping pro- 
cesses by the second order perturbation^ Each vana- 
dium ion V 3+ has d 2 configuration, and it is in a high spin 
state, 5 = 1, due to large intra-atomic interactions. Each 
vanadium site is represented by a spin state Sf — 0, ±1 
and an orbital configuration {n.i Q } 3 =1 , which is subject 
to the local constraint, ^ Q n ia = 2, imposed by the va- 
lence of the vanadium ion. Different local states are hy- 
bridized by the hopping processes, and this is described 
by an effective spin-orbital model of Kugel-Khomskii type 
on the pyrochlore lattice. It reads 
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) + (l — ^ia(ij) ) n ja(ij) , (3) 
#o-f = -KqC(1 - Si • Sj)n ia ^j)n ja ^j), (4) 

where K = t 2 a /U > 0, A = (1 - §7?)/(l - 2rf), B = 
§77/(1 - 2ry), C = (1 + rf)/(l + 2rj), and <q = J/U as de- 
fined before. ct(ij) denotes the orbital in which electron 
hopping is possible between the sites i and j. 

The two parts, H^ l _ AF and H^_ F , represent the orbital 
"antiferro" and "ferro" interactions, respectively. The 
sign of the spin exchange term in these two indicates that 
an antiferro-orbital configuration favors a ferromagnetic 
spin state, while a ferro-orbital configuration favors an 
antiferromagnetic one, as is usual for the Kugel-Khomskii 
effective model. It is noted that the orbital parts contain 
only density-density interactions as a consequence of the 
above-mentioned character of the hopping integrals. This 
means that the orbital interaction has a large anisotropy 
and the anisotropy axis depends on the bond direction of 
two neighboring sites, in contrast to the spin space with 
full rotational symmetry. 

We now give a simple argument to discuss the charac- 
ter of orbital and spin fluctuations at high temperatures 
where no symmetry breaking takes place. To discuss spin 
fluctuations, we replace orbital density operators by their 
mean value, n.i a — > (rii a ) — |, and the result is a sim- 
ple Heisenberg model on the pyrochlore lattice with the 
nearest-neighbor coupling J s — ^Kq(C — B) that is anti- 
ferromagnetic for the realistic value of 77 ~ 0.11. This is a 
strongly frustrated system and the enhancement of spin 
correlations with decreasing temperature will be quite 
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FIG. 2: Mean-field ground states of the effective orbital 
model Eq. JJJ. Configuration in a tetrahedron unit cell is 
shown. Labels, xy et al., indicate two occupied orbitals at 
each site. Solid (dash) lines denote antiferro-orbital (ferro- 
orbital) bonds, which have ferromagnetic (antiferromagnetic) 
spin couplings. 
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small. The situation in the orbital part is different due 
to its anisotropy. This time, we replace the spin opera- 
tors by their mean value, — ► (Si) = 0, and obtain the 
following effective orbital Hamiltonian, 

H ar h = K \{2A - C)n ia (ij)nj a (ij) 

)] , (5) 

but the linear terms in rii a become a constant after tak- 
ing the summation over bonds because of the local con- 
straint. Since 2A — C ~ 1 for small 77, the orbital inter- 
actions are "antiferromagnetic" , but the number of local 
orbital states is three not two, which corresponds to a 3- 
state Potts model (or equivalently a 3-state clock model) . 
More importantly, the anisotropy axis varies from bond 
to bond depending on its direction. 

We now consider the stable orbital configuration of 
izorb- As for a tetrahedron unit, there are essentially two 
different types of stable configurations shown in Fig. 2. In 
the first type (a), two ferro-orbital bonds do not touch 
each other, while they touch at one site in the second 
type (b), and these two types have the same energy, 
— (4A+2C)K . Aside from the examples shown in Fig. 2, 
there are 2 and 3 other equivalent configurations belong- 
ing to the first and second type, respectively. 

The degeneracy in energy between the two types of 
configurations is lifted if the Jahn- Teller coupling is taken 
into account. Below the higher transition temperature 
T c i, the system is compressed along the c axis. Con- 
sidering the corresponding shift of oxygen atoms along 
the c axis, the energy level of d^-orbital is pushed down 
relative to the other orbitals, and we write the level sep- 
aration as Ajt- The energy correction is then — |Ajt 
and — 5 Ajt for the two orbitals configurations shown in 
Fig. 2 (a) and (b), respectively. The energy difference 
is due to the difference in the number of occupied d xy 
orbitals. Therefore, it is natural to understand that the 
higher-temperature transition is an orbital ordering of 
the type (a). 

This orbital order introduces a spatial modulation of 
spin exchange couplings, and leads to the reduction of the 
spin frustration as a consequence. As seen from Eqs. P|)- 
a ferro-orbital bond has an antiferromagnetic spin 
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FIG. 3: (a) Spin exchange couplings in the orbital ordered 
state. Strong antiferromagnetic couplings and weak ferromag- 
netic couplings are shown by black and white bonds, respec- 
tively. J3 is the third-neighbor interaction, (b) Magnetic unit 
cell of the q = (0, 0, 2n/c) state. 



coupling, (Jia = KqC), while an antiferro-orbital bond 
has a ferromagnetic coupling (J±p — —KqB). The spatial 
pattern of spin exchange couplings is shown in Fig. 3(a). 
Since the antiferromagnetic couplings are much stronger 
than the ferromagnetic ones, Jia S> \Jif\ for the realistic 
value of 77, antiferromagnetic spin alignment is stabilized 
in each chain in the xy-planes, and the entire spin config- 
uration will be built up by stacking of antiferromagnetic 
chains. 

Interactions between the antiferromagnetic chains are 
ferromagnetic ones, Jip, but frustrated as shown in 
Fig. 3(a). The frustration is due to the structure of the 
pyrochlore lattice and also antiferromagnetic spin cor- 
relations in each chain, but the sign of interchain cou- 
plings Jif is not essential. Therefore, the relative angle 
of spins between different chains remains undetermined 
in the mean-field level approximation for the spin-orbital 
model H so , and it is determined by some other mecha- 
nisms. 

The first mechanism to consider is longer range ex- 
change interactions, in particular, the third neighbor in- 
teractions, J3. Fitting of the band calculation results 
predicted that the er-bond of third-neighbor pairs has a 
larger amplitude than for second-neighbor pairs^* and, 
more importantly, exchange couplings are frustrated be- 
tween second-neighbor pairs. The third- neighbor ex- 
change coupling has a large amplitude, and it is also 
antiferromagnetic. Therefore, this implies an order with 
q = (0, 0, 2tt/c). Two tetrahedron unit cells in two neigh- 
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J 1F /J 1A =-o.i 




FIG. 4: Zero-point energy of quantum fluctuations as a func- 
tion of 6, the angle of local staggered moments between the 
neighboring chains. 

boring xy-planes have spin configurations opposite to 
each other. However, the relative angle, 9, between the 
spins in the bottom plane of each tetrahedron and those 
for the top plane is not yet determined. 

The second mechanism is fluctuations of spins, and we 
here determine the relative angle 9 which minimizes the 
the zero-point energy of quantum fluctuations^ We use 
the standard spin wave approach starting with the mag- 
netic unit cell containing 8 sites shown in Fig. 3(b). By 
means of the Bogoliubov transformation, the energy dis- 
persion of magnons is calculated, lo^, where 7 labels the 
magnon branch, and the zero-point energy is obtained 
by Eo pt (9) = Ylkf U> ^yi where ft is the number of 
sites and h—1 in our units. Calculations have been per- 
formed for various values of Jif/Jia, and typical results 
are shown in Fig. 4. Here we set J3 = for simplicity, 



since this term is not essential for determining 9. The 
zero-point energy is always minimum at 9 = 0, 7r, which 
means that collinear order is stablest, and these two are 
equivalent, since they are related to each other by the 
mirror symmetry with respect to the x — y plane. This 
spin order agrees with the experimental result shown in 
Fig. 1. 

In this paper, we have proposed a scenario for the two 
phase transitions in the vanadium spinels AV2O4 (A=Zn, 
Mg, Cd). In our scenario, the high-temperature transi- 
tion at 50 [K] is an orbital order assisted by the Jahn- 
Teller distortion. This orbital order introduces spatial 
modulation of spin exchange couplings. The geometri- 
cal spin frustration is consequently relaxed, and this in- 
duces the low-temperature magnetic transition at 40 [K]. 
If there is no orbital order, the system remains subject 
to strong frustration, and we believe that this is the case 
for ACr 2 04 (A=Zn, Mg, Cd). In those compounds, since 
Cr 3+ ion has the (3d) 3 electron configuration, it has a 
spin S = 3/2 and no orbital degrees of freedom. Those 
compounds show only one transition at ~10[K], and we 
believe that this is driven by a coupling to lattice dis- 
tortion as proposed in Ref. [I3 • A similar transition can 
occur at a low temperature in the vanadium spinels in 
principle. However, since the orbital order and subse- 
quent magnetic transition occur at higher temperatures, 
this type of spin Jahn- Teller transition does not occur in 
the vanadium spinels in reality. 
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